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^ . Abstract 

, We present a criterion for a shock wave existence in relativis- 
tic magnetic hydrodynamics with an arbitrary (possibly non-convex) 

■ equation of state. The criterion has the form of algebraic inequality 

O , that involves equation of state of the fluid; it singles out the physical 

I solutions and it can be easily checked for any discontinuity satisfying 

' concervation laws. The method of proof uses introduction of small vis- 

Q_|. cosity into the coupled set of equations of motion of ideal relativistic 

Q I fluid with infinite conductivity and Maxwell equations. 

1 Introduction 

> ' 

^ , It is well known that not every discontinuous solution that can be formally 

^ I obtained from equations of hydrodynamics is physically allowed. The most 

well known example is prohibition of rarefaction shock waves in media with 
a convex equation of state as ones through which entropy is decreasing. But 
from the point of view of hydrodynamic equations only, such shock waves 
have the same "rights" as usual compression shock waves. This example says 
that for a shock wave to be physical some additional (non-hydrodynamic) 
restrictions must be applied. 

In case of usual hydrodynamics (without magnetic field) there are three 
famous types of additional criteria that all should hold. 

• Entropy condition. This criterion means that entropy of a fluid must 
increase as it is crossing a shock front. 
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Figure 1: Shock wave viscous profile. Dotted line: limiting case of zero 
viscosity 

• Evolutionarity conditions: Vahead > Cg, Vbehind < Cg. This criterion fol- 
lows from condition that small perturbations must be uniquely defined 



• Existence of viscous profile. This condition takes into account that 
shock wave is not a sharp step but, due to dissipative processes, some 
continuously smeared out region that we call a profile or structure (see 
Fig^. From mathematical viewpoint the shock is a generalized solu- 
tion that may be represented as a certain limit of viscous flows in zero 
viscosity limit 0. 

If all these conditions hold then one can be sure that a shock wave exists. 
In case of "normal" media these criteria are equivalent. Conditions for normal 
media are regulated by the set of Bethe- Weyl conditions f? ; and among them 
the most important is the requrement of convexity of the equation of state 
(EOS), i.e. convexity of the Poisson adiabats. 

For a normal EOS all compression shock waves exist and all rarefaction 
shock waves are prohibited. However, in case of a non-convex (anomalous) 
EOS the question of existence becomes much more complicated and the cri- 
teria 1-3 are not equivalent. Rarefaction shocks, split shocks, compression 
simple waves and complex wave configurations become possible. H. Bethe and 
H.Weyl have first studied this situation for shock waves in arbitrary fluids in 
classical hydrodynamics (see, e.g.,|2]). How the conditions 1-3 are ordered 
relatively to each other is discussed in [2]; the relativistic case is considered in 
Note that criterion 3 appears to be the most effective one and it imposes 
the most stringent requirements in case of a general EOS. 

In relativistic hydrodynamics the anomalous equation of state arises in 
super-dense media particularly with phase transitions; this has applications 
in elementary particle physics (see, e.g., jSUHl) and in relativistic astrophysics 
(super-dense matter in neutron and exotic stars, gravitational collapse and 
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supernova explosions, models of gamma ray bursts). The relativistic version 
of the convexity condition is 




> 0, T = {e + p) jr? 



(1) 



s 



where n is the baryon number density, p is the pressure; e is the energy 
density, 5* is the entropy per baryon. In a neighbourhood of phase transitions 
in super-dense matter the convexity condition can be violated. The condi- 
tions for existence of relativistic hydrodynamic shocks viscous profile in case 
of a general EOS have been investigated in OIHI. 

The situation in magnetohydrodynamics is more complicated, because 
of additional degrees of freedom due to magnetic field. In present paper 
we review the main points of investigation of existence criteria for magne- 
tohydro dynamic shocks. We introduce a viscosity term into equations of 
relativistic magnetohydrodynamics in order to study shock structure and de- 
termine conditions for a viscous profile existence. We derive these conditions 
for arbitrary EOS on basis of Landau-Lifshits viscosity tensor with one of 
viscosity coefficients put equal to zero. We show that there is a domain of 
avoidance that makes our conditions more stringent than evolutionarity con- 
ditions. Then we extend this analysis to both non-zero viscosity coefficients. 
Our consideration relies upon the results of (Zj-IHI- 

2 Basic Equations 

We consider ideal relativistic fluid with infinite conductivity in magnetic field 
with energy-momentum tensor 



where - four velocity, F^y -electromagnetic field tensor, g^" = ^^,^=diag(l,- 
1,-1,-1), h^' = -ie'^^^^F^^M^ (4- vector of magnetic field), e°^T"^ - absolutely 

anti-symmetric symbol, p* = p + ^\hf, e* = e + ^\hf, \hf = —h"ha > 
0, /i = const - magnetic permeability. The tensor Q is the sum of corre- 
sponding tensors of hydrodynamics and electrodynamics. 

Then we introduce the Landau-Lifshits viscosity tensor ^J: 

T^,u = V{u^,u + Uu,^ - U^u'^Uy^a " UyU°'u^^a) + - ^T] / 3)u"^{g - U^Uy), (3) 
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into equations of motion 



9^(r'^'^ + r'^^) = 0, (4) 

d^{nu^) = 0, (5) 

- u''h^') = . (6) 

Here Eq. (jH) describes energy-momentum conservation, Eq. (jSj) - baryon 
number conservation; Eq.(jni) follows from Maxwell equations for the electro- 
magnetic field. 

In order to obtain shock structure in the rest frame of the shock we have 
to take a one- dimensional stationary solution of (jH)-® depending upon the 
only variable x corresponding to the direction perpendicular to the shock 
front. This yields 

rplu I ^Iv 

(0)' 



u'h-' - h'u'' = H^^ u\,^h1,^ - h\,^u\,) (8) 

riM^ = n(o)M(o) (9) 

If we denote the state ahead of the shock transition by index "0" then 
relation between states on both sides of the shock (with state " 1" being the 
state behind a shock) must be as follows 

= Tlo), (10) 

~ ^tl^fel = = ^^(O)^© - ^(0)<0). (11) 



u 



n 



O^ta = ^(0)^(0) • (12) 
By definition we say that shock wave "0"^"1" has a viscous profile if 
there is a continuous solution of ©-Q having corresponding asymptotics 
for x—*oo and x—>- -oo. 
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3 The case of one viscosity coefficient 77=0. 



One can choose the reference frame such as m'^ = and h?. Then one can 
obtain 

e*u^ = Tl^^u^ (13) 

from ^ and 

= ^[H^' -u^'H''ua\ (14) 

from (jHl). 

From continuity equation 

n = U(o)no/u\ (15) 

then Eqs. fll3|) - (fTBj) allow us to express all the variables in terms of and 

Using this one can get then relationship between and 

{h\^ - H\^) {h\' - H\^) = Attu' {Tl,y - T^^nO) , (16) 

This allows us to relate three-dimensional velocity components Vi = /vP 
and V2 = /vP with each other. Note that we do not consider here the trivial 
case, when the magnetic field is zero. We have an explicit dependence Vi{v2) 
from (fT^ that leaves us the only ordinary differential equation for v} 

that describes the viscous structure of profile. Our problem is reduced to find 
a necessary condition for a regular solution of Eq. lfTTj) to exist. We assume 
also that hydrodynamic parameters in the shock structure are monotonous 
functions. If we suppose that such solution exists then we have finite dvi/ dv2 
and for monotonous structures this cannot change its sign. Then the function 
V2{vi) is well defined during the transition "0"— s>"l". In this case the right 
hand side of (fT7j) must not change its sign; otherwise the hydrodynamic 
parameters could not reach the point "1". This yields necessary condition 
we are looking for. For details see [UIH]. 

In view of we can express Vi and f 2 in terms of the specific volume 
V =l/n. 
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If introduce the function 

P{V) = {1 + (u'fr'iTl^) + i^iH-u^f - Tl^^u^}- . 

where e{V) = T^Q^Ufj_/u^ — ^\h\^ , then we obtain following criterion of ad- 
missibility of stationary shock transition: 

{V, - V,) {p{V,e{V))-p{V)) > (18) 

for all V between Vq and Vi. 

This criterion has the following consequences [7,8]: 

• If the shock satisfies the criterion, then it also satisfies the entropy 
criterion. 

• The criterion can be written in terms of the shock adiabat pu {V) as 
follows 

{Vi-Vo){ph{V)-p{V))>0 . (19) 

In case of a non-single valued shock adiabat the criterion can be for- 
mulated in terms of absence of intersections of pn {V) and p {V) every- 
where except Vo and Vi. 

• In the neighbourhoods of initial and final points we get that the pos- 
sible fast shock transitions are found to satisfy the following relations 
between the speed of the shock and the characteristic speeds at initial 
state "0" (ahead of the shock) 

VshiO) > Vf{0) > Vn{0) > Va{0) > Vsl{0) (20) 

and at the final state "1" (behind the shock). 

> "^sh^ > ^iv|2J > ^a|2J > ^sim ^21) 

Here we have introduced a new characteristic velocity that satisfies equa- 
tion 

(P + ^)Kf[K)' + l] = ^(/^')' (22) 

It is distinct from Alfven velocity that enters usual evolutionarity condi- 
tions. 
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For slow shocks we have the usual evolutionarity conditions. 

VA(0) > VshiO) > UsliO) ; Vf^> > > ^s«|2J > ^s/^H' (^3) 

So we have obtained a domain of avoidance behind the fast shock [vn, va] 
that satisfies evolutionarity conditions but nevertheless must be considered 
as non-physical one as according to our criterion there is no shock wave with 
viscous profile there. Note that the existence of f at is due to monotonicity of 
the dependence fi(f2). 

This relations for the final state appear to be more restrictive than the 
standard evolutionary criterion that does not involve Vn ■ However in the 
nonrelativistic limit V^- Va and these inequalities tend to the standard 
evolutionarity conditions. 



4 Two nonzero viscosity coefficients 



In this case when f] ^ instead of algebraic equation (jl5j) we have the second 
ordinary differential equation (for -u^). It is convenient to introduce a new 

2 



variable v = u j y 1 + {u^Y- 

After this we have the following dynamical system regarding to v} and v 



{^ + 4r^/3) — = F,{u\v) 
dv ^ , 1 , 



(24) 



where 
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1\2] 



[1 + (m1)2]3/2^,1 

Let the state parameters u^^y /i^q-,, wq, Po ahead of the shock and u^-^y h^^y ni,pi 
behind the shock satisfy the conservation laws (fTm) -()12j ) that relate hydro- 
dynamic quantities on both sides of the shock. We denote y = , yo = u^q^ 



7 




Figure 2: Fast shock. The arrows show direction of the phase trajectories 
when crossing the curves Vi and V2. Red dashed hne corresponds to viscous 
profile of the shock. 

yi = In this section we deal with the state variables u^, h^, n, p ahead 
of the shock unless otherwise stated, and we omit further the index "0" for 
these variables. 

Now we rewrite the conditions of viscous profile existence in terms of the 
right hand sides of the system and the curves Vi,V2: 
v = V,iyy. F,{y,V,iy)) = 

v = V2{yy. F2{y,V2{y)) = 

vi = Viivi) = V2{yi); vo = Vi{yo) = 1^2(1/0), Vi{y) ^ ^2(2/) Vy G (yuyo) 
The conditions can be formulated as follows 

A The function v = V2{y) is monotonous on {yi, yo). 
B For y G {yi,yo) the following inequality is valid: 

iyo-yi)Fiiy,V2iy))<0 , 

C We suppose that h^h'^ 7^0 at the point "0" (this is a technical requirement). 

D (Dl): > Uf ahead of the shock at the point "0" and un < < Uf 
behind the shock at the point "1", or (D2): ua > > Ugi ahead of 
the shock at the point "0" and Ugi > behind the shock at the point 
"1. 

The last inequalities are the famous relations between the velocities ahead 
of the shock. The first inequality (Dl) corresponds to the fast shock and the 
second one (D2) - to the slow shock. In case rj=Q inequalities (Dl) and (D2) 
follow from criterion of viscous profile existence of the section (0). 
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Figure 3: Slow shock. The arrows show direction of the phase trajectories 
when crossing the curves Vi and V2. Red dashed hne corresponds to viscous 
profile of the shock. 

Under these conditions one can show, e.g., that in case Dl the rest point 
"1" of the system (fT7|) is a saddle point. In case D2 the rest point "0" is a 
saddle point. This enables us to restore the qualitative behaviour of solutions 
to the system inside the domain of the phase plane between the curves 
Vi and V2. 

Typical situation is shown on Figs. 2 and 3. In case of the fast shock (Dl) 
the solutions leaving " 0" go out of the domain, except the only solution rep- 
resenting the separatrix of the saddle point "1"; just this solution represents 
the shock viscous structure. 

In case of the slow shock (D2) there is the only solution leaving "0" 
that tends to "1"; just this solution represents the shock viscous structure. 
The other solutions enter the domain, crossing V2 from right to left and Vi 
bottom-up. 

As the first result we have obtained following sufficiency of conditions for 
r/>0, e>0. 

Let the states "0" ahead of the shock and "1" behind the shock satisfy 
the conservation equations. If the conditions (A-D) are satisfied, then the 
MHD shock transition "O"—)-"!" has a viscous profile satisfying the equations 

More detailed proof will be published elsewhere [H]. 
If additional limitations on EOS (e.g., convexity) are absent, the criterion 
(B) dealing with EOS for the whole interval between the states "0" and "1" 



9 



is evidently more restrictive than, e.g., evolutionarity conditions [9] or any 
other conditions that involve characteristics of the fluid only at initial and 
final states. On the other hand, the evolutionarity conditions can be derived 
from the viscous profile existence [SI IH] • 

Our criteria can be applied to an arbitrary smooth EOS. However, we 
must note that the requirement for Vi{y) to be a continuous (single- valued) 
function is not trivial and can not be fulfilled in case of a certain equations of 
state (see, e.g., |21E|)- Though consideration of a viscous profile seems to be 
rather effective for investigation of shock existence and stability, this method 
can not work in case of complicated EOS (see, e.g., remarks in [3] in case of 
relativistic hydrodynamics) that require either modification of the equations 
of motion or using additional physical information about solutions. 
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